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The neutrino-driven wind from a nascent neutron star at the center of a supernova expands into 
the earlier ejecta of the explosion. Upon collision with this slower matter the wind material is 
decelerated in a wind termination shock. By means of hydrodynamic simulations in spherical 
symmetry we demonstrate that this can lead to a large increase of the wind entropy, density, 
and temperature, and to a strong deceleration of the wind expansion. The consequences of this 
phenomenon for the possible r-process nucleosynthesis in the late wind still need to be explored 
in detail. Two-dimensional models show that the wind-ejecta collision is highly anisotropic and 
could lead to a directional dependence of the nucleosynthesis even if the neutrino-driven wind 
itself is spherically symmetric. 
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1. Introduction 

The site or sites where the r-process nucleosynthesis takes place are still unclear. To be a major 
source of r-process elements, an environment must fulfill two requirements: Firstly, it has to repro- 
duce the r-process abundances pattern, and secondly, it should yield the total amount of r-process 
material in the present Galaxy. Due to this second constraint, several scenarios (e.g. neutron star 
mergers, GRB, disks) are not very likely to be the main site, although r-process elements are pos- 
sibly produced there [ffl, 0, H]. On the other hand, models for the core collapse supernova scenario, 
which was proposed already in 1957 by [Q] and ^ as the r-process site, have so far not achieved 
to produce extreme conditions needed to generate the heaviest elements observed. 

In principle, the supersonic neutrino-driven wind, which forms after the onset of the explo- 
sion, provides favorable conditions for the r-process - a high neutron abundance, short dynamical 
time scales, and high entropies. Several groups have studied this site by means of hydrodynamic 
supernova explosion simulation [Q, [7[], or by solving the stationary wind equations [Q, 

In the past years, parametric models based on the solution of the general-relativistic steady- 
state wind equations have been developed [|| ffl\. The relativistic treatment resulted in somewhat 
higher entropies than those found in Newtonian calculations [||]. However, also with this effect the 
common problem of all these studies remains that for "realistic parameters" (e.g. typical neutron 
star masses and radii) they do not find entropies sufficiently high for the r-process. However, none 
of these studies treated possible effects due to the interaction of the wind with preceding post-shock 
ejecta consistently. This interaction leads to the formation of a wind termination shock (or reverse 



shock) []10[ |llj], which can change the neutrino wind properties considerably. 

We have performed one- and two-dimensional hydrodynamic simulations that start a few mil- 
liseconds after bounce and follow the evolution for several seconds after the onset of the explosion. 
This allows us to study systematically the influence of this reverse shock, which turns out to in- 
crease the wind entropies significantly. Here we describe briefly the numerics we are using and 
some of our first results. 



2. Numerical setup 



Our computational approach is described in detail in [12]. The initial data at ~ 10ms after 
bounce are provided by Boltzmann simulations. We replace the inner part of the neutron star 
(p > 10 g/cm 3 ) by an inner boundary, which allows us to avoid strong time step limitations and 
provides a degree of freedom for systematic variations. The contraction of the inner boundary can 
be parametrized by a time scale and a final neutron star radius. This can be justified by the fact that 
the neutron star evolution depends on the still incompletely known dense matter equation of state. 
The core neutrino luminosity is put in "by hand" at the inner boundary (below the neutrinosphere). 
Its value is chosen such that the explosion energy reaches typical values known from observations. 
We use a simplified neutrino transport treatment (a grey, characteristics-based scheme), which 
reproduces the results of Boltzmann transport simulations qualitatively. However, quantities like 
the electron fraction should be taken with caution. Our hydrodynamics code is Newtonian, but 



we account for relativistic effects by using corrections in the gravitational potential [13]. This 



approximation has been shown to yield results very similar to the full general relativistic treatment. 
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3. Results 

After the explosion has been launched successfully, the supernova shock propagates outwards 
and the density around the neutron star decreases. Ongoing neutrino-energy deposition leads to the 
formation of a neutrino-driven wind that emerges from the neutron star surface. The matter in the 
wind is accelerated to supersonic velocities, hits the slower moving preceding ejecta and is strongly 
decelerated in a wind termination shock (see Fig.[|). 
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100 




Figure 1: Mass shell and entropy (color) plot for a one-dimensional model (based on a 15M© progenitor) 
that attains an explosion energy of ~ 1 .2 bethe. The green line corresponds to the neutron star which shrinks 
to a radius of ~ 10 km. The red line is the shock radius. The dashed lines are mass shell trajectories (the 
baryonic mass below a mass shell is given in M Q ). In the wind phase the mass shells move to larger radii 
in a very short time corresponding to high wind velocities. When the wind velocity becomes supersonic in 
the frame of the slower-moving material outside (i.e. in the region where the mass shell lines expand more 
slowly), a reverse shock forms. This leads to a jump in the entropy (at ~ 5 • 10 8 cm). 

Figure || shows characteristic quantities from a one-dimensional simulation as functions of 
radius for different times. In the wind phase, the velocity increases, becomes supersonic, and tends 
to an asymptotic value. Also the entropy reaches a constant maximum value where the neutrino 
energy deposition rate (q in Fig.Q) becomes negligible. The quantities characteristic for the wind 
(velocity, entropy, ...) are qualitatively and, for the same values for the neutrino heating rate, even 
quantitatively in agreement with general relativistic stationary wind solutions (see e.g. The 



wind entropy increases with time due to the decrease of the neutrino luminosity (see [[14)] for an 
analytic discussion of how the entropy depends on the neutrino properties). The highest wind 
entropy obtained in our simulations is not as high as those found in previous wind studies. The 
reason is that the neutrino luminosities even at the end of our simulations are higher than the lowest 
values considered for stationary wind solutions. This explains also the larger mass flux in the wind 
(M lO~ 4 M /s) in our case. Another nucleosynthesis-relevant quantity, the electron fraction 
Y e , evolves from proton richness (Y e > 0.5) just after the onset of the explosion (as found also in 
Boltzmann simulations with artificial explosions [^, [16|]) to neutron richness (Y e < 0.5) at t > 3s 
after bounce. In the simulation presented here we have short expansion time scales which should 
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prevent a large fraction of alpha particles from recombination to seed nuclei. Consequently the 
situation favors many free neutrons to remain unbound and thus a high neutron-to-seed ratio after 
the a -rich freeze-out phase. 




Figure 2: Radial profiles at different times after bounce for the neutrino energy deposition rate (q), density, 
entropy (left column), velocity, temperature and 5^ (right column). The effect of the reverse shock is clearly 
visible as a jump of the wind entropy by a factor of two. 



In Figure ^| one can also see how the supersonic outflow hits the slower preceding supernova 
ejecta, and is decelerated in the wind termination shock, which leads to a sudden increase of den- 
sity, temperature, and entropy. From the Rankine-Hugoniot conditions it is straight forward to 
demonstrate that the jump in the entropy increases with the square of the wind velocity. Tempera- 
ture and density decay with time much more slowly after the wind matter was decelerated by the 
reverse shock. The conditions established during that phase (T9 ~ 0.5 — 1, p ~ 100 — 10 4 g/cm 3 ) 
suggest possibly important consequences for the nucleosynthesis in the supernova ejecta. Detailed 
nucleosynthesis calculations, however, are needed to explore the exact effects. 

We have performed simulations of the kind described here, with varied conditions at the inner 
boundary (different neutron star contractions and neutrino luminosities) for different explosion 
energies and several progenitor stars. In all cases the wind termination shock developed and had 
qualitatively similar consequences. Interesting quantitative differences will be analyzed in a paper 
in preparation. 

The reverse shock forms also in our two-dimensional simulations (see Fig. || for an example). 
Interestingly, the conditions are strongly dependent on the direction and suggest a very anisotropic 
distribution of the nucleosynthetic products of the neutrino-driven winds. This also raises the 
question how robust the condition can be in this environment for r-process nucleosynthesis. 
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Figure 3: Two dimensional simulations. Entropy and velocity at 1.2 seconds after core bounce. Note the 
highly aspherical distribution of the ejecta. The effect of the reverse shock on the entropy leads to similar 
maximum entropies as in one dimension. 



4. Conclusions 

Our work demonstrates that over a wide variation of conditions (neutron star parameters, pro- 
genitors, spherical symmetry or multi-dimensional environment) the wind termination shock may 
have non-negligible influence on the nucleosynthesis conditions in the neutrino-driven winds. It 
does not only increase the wind entropy by up to a factor of two, but also changes the time evo- 
lution of density and temperature in an interesting way. Nucleosynthesis calculations are needed 
to study the exact consequences of this so far incompletely explored feature of supernova ejecta 
dynamics. 
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